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1 Multilevel determinantal ensembles 



1.1 Joint probability distributions 

The following common setting underlies various multilevel determinantal ensembles, in- 
cluding chains of random matrices [2j[3l[8], and finite determinantal point processes such 
as Dyson processes [10] and polynuclear growth [6]. Let {(Tj,dfij)}j=i... m be a set of 
measure spaces and {Hj := L 2 (Tj, rf//j)}j=i...m the Hilbert spaces of square integrable 
functions on them. Suppose we are given m — 1 functions {wj+i, J } ;? =i...m-i on the product 
spaces Tj + i x Fj, such that the corresponding integral operators 

v>i + iAf)(* U+1) ) ■= I w J+1 j{x<* +i >,x<»)f(x®)dto(x®), (1.1) 
together with their transposes 

w* +ld (f)(x^) := I w j+ld (xV +1 \x®)f(xU + Q)d^xV +1 '>), (1.2) 
Jr j+1 

are well defined injective maps on a sequence of dense subspaces H± C Hi , H 2 = 
W2i(Hx) C H 2 , . • • , i? m = w mjm -i(H m ..i) C # m , as are their composites: 

w kj := w k ^io . . .o Wj+1Jl m>k>j>l. (1.3) 

Let ifi C -f/i be an N dimensional subspace with basis {ip^} a =i,...N, and {Hj : = 
Wji(Hi) C Hj}j =lj jm , the corresponding subspaces obtained by applying the operators 
{wji} to ifi, with bases 

ty® := Wjl (^)} a=1 ,... N . (1.4) 
Now choose the basis {<j>a} a =i,...N f° r dual to {V'i }a=i,..JV> and let 

#:=<(0i m) ), a = l,...iV (1.5) 

be the corresponding dual bases for the Hj's. Thus 

/ ^\x {j) )^\x^)d^{x^) = 5 ab . (1.6) 
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Assume that 

m— 1 m N 

det(^(x 5 w ))det(^(xl m) )) J] detK. +1)i (^ +1 ) ! 4 J ' ) ))nn^^' ) ) ( LT ) 

j=l j=l a=l 

is a positive Borel measure on n^li^i)^ It follows from the duality relations (11.61) that 
this is a normalized probability measure 

n / nc^s )) det (^ i) (4 i) ))det(0i m )(xi m) )) 

j=i V/rf 0=1 / 

m— 1 

xHdetK^fxfVfVl (1.8) 
i=i 

so that 

m—l 

P N > m (x^) := det(^ 1) (4 1) ))det(0i w) (4 m) )) U det(w s+ltj (x^ +1 \ x®)) (1.9) 

j'=i 

may be interpreted as a joint probability density for points in OJLi (Tj) 7V - 

A typical example of such probability measures arises in the theory of random matrices. 
We start with a chain of N x N Hermitian matrices {Mj 6 M. N }j = i m with nearest 
neighbor couplings of exponential type e tr ( M i M i+i) ; anc l each site in the chain contributing 
a conjugation invariant factor e~ tT ^ Vi ^ Mi '> to the joint probability density (where the V/s 
are e.g. polynomials in the matrices), and the product gives a normalizable positive 
measure. Integrating over the "angular variables", using the Harish-Chandra-Itzykson- 
Zuber (HCIZ) identity [5], leads to a reduced joint probability measure of the form ( 11.91) 
on the space of eigenvalues, where the T/s are all taken as the real line, 

w J+lJ (x^ +1 \x^) := ^H^'HVi+i^B (L10 ) 

V>iV 1} ) := Pa -i(x (1) )e-^ m \ 0i m) (^ (m) ) := Sa^x^e-^^ (1.11) 

and {p a (x^), s a (x^)} a= o,i,... are a sequence of pairs of polynomials of degrees a satisfying 
the biorthogonality relations 

m f n \ m—l 

II y j vN d ^ ] )j ^K^)4 n \x^) J! w j+ld (xV +1 \x®)) = 5 ab . (1.12) 

If all the iys are identified as the same space T, (\1.9\i may also be interpreted as a 
measure on the space of m-step paths of iV-tuples of points in T, defining a point process, 
in which j is viewed as a discrete time parameter. Examples of this kind lead to Dyson 
processes [10], describing diffusion of eigenvalues, and Polynuclear Growth [6]. 



3 



1.2 Point correlators and gap probabilities 

What characterizes such multilevel determinantal ensembles is that multi-point correlators 
(marginal distributions), gap probabilities and expectation values may all be expressed 
as determinants, either finite, or determinants of Fredholm integral operators, in terms of 
a single m x m matrix kernel function Kij(x^ l \x^')i ) j ) =i l ... m defined as follows: 

Kijix^x^) := K tj {x^,x^) - Wij (x^,x^), (1.13) 

where 

N 

K i:i (x®,x®) :=Y,^a ) ( x(i) ) ( Pa ) ( x{3) )- (1-14) 

a=l 

wyix®,!®) :=0 if %<] (1.15) 

and Wij(x^,x^) is defined as in (II. 3ft for % > j. In terms of this, the joint probability 
density (jl.9p may be expressed as 

P N >™{xf)=det{k ij {xf,x ( i ) )) (1.16) 

where Kij(xa\ x^ ) is viewed as the ((z, a), (j, 6)) element of a matrix of dimension iVm x 
Afm, labelled by pairs of double indices 1 < i,j < m, 1 < a, b < N. By integrating eq. 
(I1.16p over some of the variables while setting the rest equal to fixed values, it follows 
[2] that the correlation function giving the probability density for finding kj elements in 
Tj at the points {x± , . . . x^ } for j — 1 . . . m is similarly given, within a combinatorial 
factor, by the Y^T=i x Y^j=\ determinant 

/&, ({•''," %h ,..,,/) : det(J^(4« af)) , .. , . (1.17) 

Alternatively, choosing a measurable subset J,- C Tj of each Tj, the probabability E(0, J) 
of finding no points within the set 

J := Ji X J 2 X . . . J TO (1.18) 

is given, within a normalization constant, by the Fredholm determinant ([8], [TO] ) 

^' m (0,J) = C iVm det(J-/Co Xj ), (1.19) 
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where K, is the mxm matrix integral operator with kernel Ky defined in (11.131) and xj is 
the direct sum of the operators of multiplication by the characteristic function {xjAj=i,...m 
on each of the factors in the direct sum 

H := @f =1 Hj. (1.20) 

More refined statistics, giving the probability of finding any specified number of points 
within disjoint subintervals {Jji}i=i imj of the Jj's may be similarly computed by replacing 
the characteristic functions \Jj by weighted ones ^a=\ z jiXJ t an d evaluating the coeffi- 
cients of the monomials Y\T=i nl=i z *jf ■ 

A determinantal expression similar to (I1.17P may also be found for the so-called 
"Janossy densities" [H [8] , which give the probabilities of finding elements at the given 
set of points {x±, . . . x^} within the subsets (Ji, J 2 , . . . J m ), while all others are on the 
outside, normalized to the probability of all points being on the outside. In that case, 
the integral kernel Kij(xa\x^) appearing in (11.171) needs simply to be replaced by the 
corresponding kernel Rf?(xa\ ) of the Fredholm resolvent operator: 

^^(l-^o^, (1.21) 

where 

/C^:=/Co Xj . (1.22) 

These results may all be obtained as consequences of a single determinantal identity 
for multilevel ensembles, which can be expressed in a purely algebraic form. In the next 
section this algebraic identity will be stated, and a very simple proof given. All the 
aboveexpressions for correlators and gap probabilities may be deduced from it, as can 
similar expressions for arbitrary multilevel determinantal ensembles. 

This result is not at all new; it was derived e.g. in refs. [2], [3], [8], [10] in a number of 
particular cases, using a variety of different methods. The point of giving a new proof of 
the underlying identity is just to unify and simplify the argument, showing its universal 
applicability. By suitably particularizing to the various cases of interest, all previous 
results of this type follow, and determinantal expressions may be similarly derived for 
further cases that have not previously been studied. Some examples, consisting of chains 
of Lie algebra valued matrices, including those having various nearest neighbour couplings 
that admit reductions similar to that of Harish-Chandra-Itzykson-Zuber, are given in 
section 3. 

Following Tracy and Widom [9j [10], the relevant integrals are first obtained in a 
uniform way by introducing a set of m measurable functions or distributions {pi, . . . , p m } 
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on {TV, ... , T m } and considering the integral obained by replacing each dp,j by dp,j{l — pj). 



N 



P ZZ Pm ■= II / II (^i(^)(l - Pj(*?)) det(^ 1 )(xi 1) ))det(^(xi m) )) 

171 — 1 

x l[det(w j+1J (x^ +1 \x^)) (1.23) 
j'=i 

We then apply a multi-level version of the Andreief identity PQ to express this integral as 
an N x N determinant. 

P pi ,..., Pm := (iV!) m detG, (1.24) 

where 

m / » \ m—1 

(1.25) 

Like the one-level version, this follows easily from the invariance of the integrand under 
permutations in the various integration variables at each level. 

By choosing the factors {pj} in different ways, Tracy and Widom [9], [10] showed that 
point correlation functions, gap probabilities and Janossy densities may all be obtained 
from the formula fl 1.23ft . For example, choosing the pj's as: 

kj 

pj := - ZjiS x u) (1-26) 



i=i 



the coefficient of the monomial factor YlJLi Yl^Li z ji i n Ppi,...,p m gives the point correlator 
Pk x> ...k m {{X\ J \ . . . Xkhj=i..jm)- Choosing instead 



i=i 



gives the generating function for probabilities of finding the various numbers of points 
within disjoint subintervals Jji C Jj. Choosing 



Pr-=Xj s -^2zji6 xt (1.28) 



i=i 
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and evaluating the coefficient of the monomial factor YlJLi Yli=i z ji i n Ppi,-,p m gi ves the 
corresponding Janossy densities. 

In order to apply the abstract form of the determinantal identity derived in the fol- 
lowing section, it is convenient to interpret the matrix elements of G as scalar products. 
To express integrals involving the functions ipa and <f$ as scalar products, these will be 
represented as vectors or linear forms, using the bra and ket notation: < ipa^\, < </>i"^|, 
> and \cj) { a j) >. Thus, G ab may equivalently be written as: 

G ab =< ^1(1 - Pl )w* 2!l (l - P2 ) ... (1 - p m _ 1 X +lim (l - Pm )\<P { b m) > (1.29) 

where 

Pj : //,— //; (1.30) 

is the operator of multiplication by pj, which is symmetric (pj = p*) with respect to the 
scalar product 

<v ij) \u {j) >f= [ v^(x^)u { - j) {x { - j) )dix j {x^) (1.31) 

on Hj. 



2 Multilevel determinantal identity 

To provide a uniform algebraic setting for all the previous examples, we assume that we 
have m Euclidean inner product spaces {Hi, . . . , H m }, and denote their direct sum by H 
as in (jl.20p . We also assume that there are m — 1 injective linear maps: 

w j+i,j '■ H i^ H 3+ii j = 1, . . . , m — 1 (2.1) 

with their duals 

w* +hj : H j+1 ^Hj (2.2) 
as well as m endomorphisms of the spaces Hj 

Pj : Hj^Hj (2.3) 

that are symmetric pj = p*. The map p : H^H is defined to be the diagonal one 
leaving each subspace Hj invariant and equal to pj when restricted to Hj. We also define 
composite maps 

w kj : Hj^>H k , m > k > j > 1 (2.4) 
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as in (11.31) . with Wkj '■— for k < j. Together, these define an operator 

W : H^H 



(2.5) 



that is lower triangular with respect to the decomposition (11.201) and can be written in 
matrix form as 



(2.6) 
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This may be expressed as 
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We also denote the corresponding 
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(2.7) 



(2i 



(2.9) 



(2.10) 



As previously, let Hi C Hi be an dimensional subspace with basis {^o }a=i,..jVj 
and {Hj := Wji(Hi) C Hj}j = i t _ >m , the corresponding subspaces obtained by applying the 



8 



operators {wji} to Hi, with bases {ipa^ 1 := w ji(ipa)} o=i,..jv) • Let {0i" 1 ' ) }a=i,...Af be the 



basis for H m dual to {-0i m) }a=i,...v, and {0 a := ^m 7 (0i m) )}a=i,...iv)} the corresponding 



dual bases for the -ff/s, so that 

<ij U)\ (f) U) >=6abi J = i,..., m . (2.11) 

For a = 1, ... N, we introduce the m-component row and column vectors 

/ > \ / > \ 

*a := ; , $a := ; (2.12) 

Vi4 m) >/ Vi0i m) >/ 

^ := (< V>«|, • • • , < ^i m) |), € := (< • • • , < 0i m) |), (2.13) 
which may be viewed as elements of H and H*: 

m m 

*.=X>N4°>, *. = X>w°> 

i=l i=i 

*2 , = E < ^Vi, C = E<^V„ (2.i4) 

3=1 j=l 

where 

', 'Hi r-^ ©7=i «i *i : ®l=iH3^% (2-15) 
are the standard injection and projection maps. These may also be expressed as 

tf a = (J-A/TVI^ >= (Z + W)^ > (2.16) 

$ Q = (X - AH~ V|0 a m) >= (1 + W) v|0i m) > (2.17) 

*I =< ^Vi^-AT*)- 1 =< ^kx(J + W) (2.18) 

C =< ^VJI-ATT 1 =< 0i m V m (X + W*). (2.19) 

The duality relations (11.61) may equivalently be expressed as 

< ^liT-N*)' 1 ^ >= S ab . (2.20) 
In terms of the elements \l/ a , <J> a G i7, define two maps 

:C N ^H, $ : C^if 

* :e a ^ $ : e a i-> $ a (2.21) 
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and their transposes 



^ T : H^C N , 



\V {1) > \ N m 



\ v (m) > J a=l j=l 



<£> T : H^C N , 

( \V W >\ N m 

$T; : hEE<i°V j) >e a) (2.23) 

V |w (m) > / a=i -? =i 

where {e a } a =i,...,N is the standard basis for for C^. 
Denote by 

/C:=^o$ T (2.24) 

the linear operator /C : H^H which is the composite map. This may be expressed in 
m x m matrix form as 

(A^ii • ■ • K>lrn \ 
: ••. : (2.25) 
^ml ' ' ' fc"mm / 

where 

N 
a=l 

We also define the linear operator K, : H^H as 

JC:=K-W (2.27) 

and let X : TL— >TC denote the identity map. Finally define, as in (11.291) the N x N matrix 
G whose components G a b are the scalar products 

G ab :=< ^1(1 - p^w^l - ^ ... (1 - p OT -i)< +lim (l - Pm )\<P { r ] > ■ (2-28) 
The main result of this section is the following identity 
Proposition 2.1 

det(G) = det(J-/Cop). (2.29) 
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Proof: We begin by denning m — 1 modified linear maps 

Wj+iJ : Hj^H j+1 , j = 1, . . . , m - 1 
by composition with (1 — Pj) : Hj^Hj 

w j+1J := w j+hj o (1 - Pi ) 
We also define, analogously to Af and W, the operators Af and W by 



W := Afo (I-Af)- 1 



where 

w, 

The dual operators are 
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Note now that 



((X + W*) o (X - p)) =71,(1 + W*) o (X - p) Lm = w* 0...0 w* o (1 - p m ) 

\ / lm 

= (1 - px)w* 2l (l - p 2 ) • • • (1 - p m -x)w* m+1 (1 - p m ), (2.37) 



and hence 



=< tfV^X - (X - P) o A/"*)- 1 o (X - P )L m \<P { r ] > • (2.38) 
From fl2T20|) it follows that 

Gafc = 5 ab - 7 a b, (2.39) 

where 

lab :=-< # K ((X - (X - p) o AT*)- 1 o (X - p) - (X - A/"*)- 1 ) ^|0i m) > 
=< V^Vi ((2 - (X - p) o AT*)- 1 o p o (X - AT*)- 1 ) U0i m) > 

= ^ ((X - A/"*) o (X - (X - p) o AT*)- 1 o p) (2.40) 

(where (12.171) . (12.181) have been used in the last line). Now applying the well-known 
identity 

det(l - A o B) = det(l - B o A) (2.41) 



to the operators 



5:=po$: C N ^H 

A:=V T o ((1 - AT*) o (X - (X - p) o AT*)- 1 ) : H^C N , (2.42) 



we obtain 



det(G) = det ((X - p o $ o ^ T (X - A/"*) o (X - (X - p) o A/"*)" 1 ) 
= det (X - (X - AT o (X - p))- 1 o (X - AT) o K o p) 
= det ((X - AT o (X - p) - (X - A/") o /C o p) 
= det ((X - (X - A/V 1 o AT o p) - /C o p) 

= det (X - (/C - W) o p) = det (l-JCopY (2.43) 
where in the third and fourth lines we have used the fact that 

det (X - M o (X - p)) = det (X - M) = 1. (2.44) 
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3 Examples 

3.1 Classical Lie algebra chains 

As furthert applications of the above identity, we consider coupled chains of m random 
matrices having values in the Lie algebras of the compact forms of the classical Lie groups 
[7j, with Itzykson-Zuber nearest neighbour exponential couplings. The partition function 
is of the form 

777, / p \ 771 771—1 

z n = n ( / dA i ) n e~ tTVj(Aj) n (3.1) 
3=1 \j / 3=1 3=1 

where {Aj e 0}j=i,... m for any of the classical Lie algebras g = u(iV), o(2N), o(2N + 1) or 
5p(N) and dAj denotes Lebesgue measure on each of these. Using the Harish-Chandra 
identity to integrate out the angular variables reduces these to integrals over the Cartan 
subalgebras. For the case u(N), identified as Hermitian matrix chains, this gives the usual 
Itzykson-Zuber reduced integral [51 [2] 



m / „ N \ 771—1 

oc J] / ]Jd^e- v ^ A(x«, . . .,x$)A(x<T\ . . .,*>>) ]J det 



3=1 

(3.2) 

where 

A(xi, . . . ,x N ) := det(x^ _a )|i< ai6 < iV = ] [ (x a - x b ) (3.3) 

l<a<ft<JV 

is the Vandermonde determinant. For the other classical Lie algebras it gives (see [7] for 
the case m = 2) 

o(2N): 

m / N \ 

Z ^ « II ( / RJV n^e^^^j A((x«) 2 , . . . , (*«) 2 )A((a^) 2 , . . . , (z^) 2 ) 

771—1 

x Y[ det(2cosh(4 i) 4 i) )) ( 3 - 4 ) 
i=i 

o(2iV + 1) or sp(iV): 

m / „ N \ / N \ 
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N \ m— 1 



A((4 m) ) 2 , . . . , (4 m) ) 2 ) J] det(2sinh(*^)), (3.5) 

\a=l J j=l 

where the potentials {Vj}j=i m must be chosen as even functions of their arguments. 

Replacing the monomial elements the Vandermonde deter- 

minants in (13.21) by (ip a (x^), (p a (x^)j, where 



:= Pa-i{x^)e-^ W \ 
4> a (x (m) ) ■= s a ^{x^)e-^ xim) \ a = l,...N (3.6) 

and p a (x^), s a (x^ m ') are polynomials of degree a G N satisfying the biorthogonality 
conditions 

m / » \ m— 1 

we obtain the joint probability density (jl.9p with (1^-, cfyty) = (R, dx^),j = 1, . . . , m and 

w ]+1 ^ +1 \x^):=e xU)xU+1) . (3.8) 

Thus, applying Proposition 12.11 to this case, and choosing the functions pj to be the 
characteristic functions xjj °f unions of subintervals Jj = \J t Jji C R gives the gap 
probability 

P"C, XJm ■= E N ' m (0, J) = C Nm det(l - £ o xj) (3.9) 
where the mxm matrix Fredholm operator JC is defined as in 

dm, (mn with (4°, 4°) 

defined in (OL (Oj) . fTL5|) . 

The other two cases are nearly identical, the only difference being that, for o(2N), the 
potentials Vj(x^') should be taken as even functions, and the biorthogonal polynomials 
(p a -i(x^)s a -i(x^)) in (13.61) replaced by even polynomials of even degrees 

MX (1) ) ■= P2 { a-l)(xU)e-^* W \ 

: = s 2 (a-i)(x {m) )e-^ xW) , a = 1, . . . N (3.10) 



'a 



X 



satisfying the biorthogonality conditions 



'''//■ \ 771 — 1 

II yj N dx {j) e- V > (xU)) J P2a{x [1) )s 2b {x^) n 2cosh(^ +1 ^)) = S ab , (3.11) 
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The functions Wj+ij entering in (11. 3p . (II .4p . (I1.5P are 

w J+ltj {x {j+1 \x U) ) := 2cosh(x (j+1) x {j) ). (3.12) 

For o(2A^ + 1) and sp(N) the potentials Vj(x^') must again be taken as even func- 
tions, the biorthogonal polynomials (p a -i(x^) s a -i(x^)) in (13.61) are replaced by odd 
polynomials of odd degrees 

^(x (1) ):=p 2cI -i(x«)e-^ 1 ^ (1) ), 

a (x (m) ) := s 2a . 1 (x^)e'^ xW \ a = 1, ... N (3.13) 
satisfying the biorthogonality conditions 

m / p \ m—1 

II y j vN (lrJ ' ( r "' j P2a-i(z (1) K&-i(z (m) ) J] 2sinh(^ +1 )^) = <U (3.14) 

and the functions Wj+ij entering in (11.31) . (11.41) . (11.51) are 

w j+1 j{x {j+1) ,x {j) ) :=2smh{x {j+1) x ij) ). (3.15) 

Remark 3.1 In view of the fact that the joint probability densities for the cases o(2N), o(2N+ 
1) and sp(N) are all even functions of their arguments, if averages only over even functions are 
considered, as would be the case, for example, in computing expectation values of conjugation 
invariant functions of the matrices, the coupling terms 2 cosh(x^- ) j;^ +1 ^) and 2 sinn^^x^ 1 )) 
appearing may equivalently by replaced by the usual exponential couplings e x x 

3.2 Nonexponential couplings 

As explained in [I], Appendix A, another generalization of the Hermitian matrix chain 
can be made by choosing a set of functions fj{z) having a Taylor series expansion about 
the origin of the form 

oo 

fj{x) = 1 + ■■■r j (k)x k , j = l,...,m-l. (3.16) 

k=l 

The exponential coupling term e tr 3 in (13. ip may be replaced by the more general set 
of conjugation invariant functions 

TjiAjAj+i) := Yl [l[r j (N + k-i))d x , N s x (A j A j+1 ), (3.17) 

e(x)<N 
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where s\ is the Schur function corresponding to partition A = (Ai > • • • > X^x) > 0) of 
length £(X), viewed as a character, and d\,N is the dimension of the irreducible tensor 
representation of (5l(N) corresponding to A. The joint probability density for the eigen- 
values on the chain is then again of the form (11.91) with the functions (^(a^ 1 '), 4> a (x^)) 
defined as in (11.111) satisfying the biorthogonailty relations ( 11.121) . but with the functions 
Wj + ij defined as 

w j+1J (x^ +1 \x {j) ) := (3.18) 

An example of such a coupling is obtained by choosing a set of m — 1 pairs of constants 
{(ai, z t ), . . . , (a m _i, z m -i)} and setting 

f j (x):=(l-z j x) N - a *- 1 , j = l,...,m-l. (3.19) 

Hence 

w j+h] (x il+1 \ x®) = (1 - Z] x^x {j+l) ) N - a >-\ (3.20) 

which corresponds to 

, . za (a,j — N + i) , 

rj(i) = ■ J ^ J —. (3-21) 

This choice gives 

r^AjAj^) = det(I - ZjAjAj+x)-** (3.22) 
and hence the partition function (13.11) for the coupled random matrix chain is 

Z n = II ( / dA i) U e ' tlVjiAj) II det(l-z 3 A 3 A J+1 )- a >. (3.23) 
j= i V / j= i j= i 

All the above formulae for multi-point correlators, gap probabilities and Janossy densi- 
ties are applicable to these cases, simply by adapting the general fomulae to these choices 
for the coupling functions Wj + ij. 
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